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Relaxed atomic geometries and chemisorption energies have been calculated for the dissociative 
adsorption of molecular hydrogen on vicinal Si(OOl) surfaces. We employ density-functional theory, 
together with a pseudopotential for Si, and apply the generalized gradient approximation by Perdew 
and Wang to the exchange-correlation functional. We find the double-atomic-height rebonded Db 
step, which is known to be stable on the clean surface, to remain stable on partially hydrogen- 
covered surfaces. The H atoms preferentially bind to the Si atoms at the rebonded step edge, 
with a chemisorption energy difference with respect to the terrace sites of 0.1 eV. A surface with 
rebonded single atomic height Sa and Sb steps gives very similar results. The interaction between 
H-Si-Si-H mono-hydride units is shown to be unimportant for the calculation of the step-edge 
hydrogen-occupation. Our results confirm the interpretation and results of the recent H2 adsorption 
experiments on vicinal Si surfaces by Raschke and Hofer described in the preceding paper. 

PACS numbers; 82.65.My, 68.45.Da 
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I. INTRODUCTION 

The interaction of hydrogen with silicon surfaces has 
become an intensively studied matter. There are impor- 
tant applications in semiconductor technology, such as 
the passivation of surfaces, etching, and chemical vapor 
deposition growth. The attachment of the deposited 
Si atoms at surface steps is an essential aspect of epitaxial 
growth. U^l^ From the difference in hydrogen chemisorp- 
tion energy between the adsorption sites on the terraces 
and at the steps of a vicinal Si(OOl) surface the equi- 
librium hydrogen occupation of the various surface sites 
can be derived. If hydrogen atoms preferentially bind 
to the step edge, already small hydrogen coverages are 
sufficient to saturate the Si dangling bonds at the step 
edge, and thereby to affect the probability of Si atoms 
to become attached to the step. In this way the kinetics 
of Si epitaxial growth can be affected by the presence of 
hydrogen. Of course there are also other mechanisms, 
like the hydrogen-induced change of the Si surface dif- 
fusion coefficient. [Q| Furthermore, within the framework 
of thermodynamic equilibrium theory, the step energies 
govern the surface morphology on a mesoscopic length- 
scale. As step energies depend on hydrogen cov- 
erage, the adsorption of hydrogen on the surface could 
affect step roughness or even surface morphology. |]l2| 

Last but not least, the dissociative adsorption and as- 
sociative desorption of molecular hydrogen on the Si(OOl) 
surface has attracted a great attention, in particular be- 
cause adsorption and desorption experiments have led to 
apparently contradictory results with respect to the ad- 
sorption energy barrier. On the one hand, the observed 
small sticking coefficient suggests a large adsorption en- 
ergy barrier, while, on the other hand, the kinetic en- 
ergy distribution of the desorbing hydrogen molecules is 



nearly therm al, sug gesting only a small adsorption en- 
ergy barrier. [y_3 14 1 Despite intensive research there are 
still competing explanations. In the dynamical model by 
Brenig et al. the desorption proceeds from two hy- 

drogen atoms bound to the two Si atoms of a single sur- 
face dimer. This hydrogen pairing is corroborated by the 
observed first-order desorption kinetics of the hydrogen, 
and the deviations towards second order at very small 
hydrogen coverage. [pT| Brenig's model contains a bar- 
rier; the small kinetic energy of the desorbing particles is 
ascribed to an efficient energy transfer into the Si surface 
phonon degrees of freedom during the desorption process. 
However, several quantum-chemical cluster-calculations 
[p"8[-pl| have arrived at distinctly larger barriers than 
density-functional slab-calculations, p^-p^ and these 
large desorption energy barriers appear to be at variance 
with the desorption dynamics sketched above. Thus the 
authors of these calculations favor models, in which the 
adsorption on and desorption from surface defects play 
a major role. [p5|-p8[ On the other hand, Pai and Doren 
pof obtained barriers from density-functional cluster cal- 
culations that are consistent with the pre-pairing model. 

To investigate the role of steps and defects experimen- 
tally, Raschke and Hofer 1 30 1 have studied the adsorption 
of H2 on vicinal Si(OOl) surfaces. For studied miscut an- 
gles larger than 2° double atomic height Db steps pre- 
vail on the clean surface. [ p^j3l| Raschke and Hofer infer 
from their experimental results that H2 molecules pref- 
erentially adsorb at the step sites. Most importantly, 
they can distinguish between the contribution of terraces 
and steps to the H2-sticking coefficient. Contrary to the 
large barrier towards H2 adsorption on the flat surface 
and on the terraces, there appears to be only a rather 
minor (^0.09 eV) adsorption energy barrier along the 
reaction path that leads to dissociative adsorption at the 
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step edge. ^ The measured hydrogen step-saturation 
coverage of about one H atom per (1x1) surface lattice 
constant a along the step is consistent with a model in 
which the hydrogen atoms bind to the edge of the dou- 
ble atomic height steps in a local mono-hydride like ge- 
ometry. Furthermore, at elevated temperatures Raschke 
and Hofcr observe the diffusion of the adsorbed hydrogen 
atoms from the step onto the terrace. From the measured 
equilibrium hydrogen coverage on the step edge and on 
the terrace they deduced a difference between the local 
chemisorption energies of roughly 0.15 - 0.3 eV by which 
the hydrogen atoms bind more strongly to the step edge 
Si atoms. [|| 

In view of the great importance of H2/Si(001) as a 
model system to understand the dynamics of adsorption 
and desorption, and in view of the conflicting results of 
quantum-chemical and density-functional theory (DFT) 
based computations of the adsorption energy barrier of 
the flat surface, it appears to be quite desirable to calcu- 
late the H2 chemisorption energies and adsorption energy 
barriers for vicinal surfaces. A successful comparison of 
DFT results to the new experimental data would lend 
support to both the interpretation of the experiment and 
the credibility of the generalized gradient approximation 
to the exchange-correlation functional used in the DFT 
calculations. 

The purpose of this work is to provide a comprehen- 
sive overview over H2 chemisorption on Si(OOl) vicinal 
surfaces: Relaxed geometries and chemisorption energies 
will be presented for various adsorption sites and step 
topologies. The adsorption energy barriers on the terrace 
and at the step edge will be discussed in a forthcoming 
paper. |32[ 

The organization of the paper is as follows: First, 
the stable configurations of the clean steps (whic h are 



a function of miscut angle) are reviewed in Section [II A 



These structures have already been thoroughly inves- 
tigated both experimentally p^-|3^ and theoretically. 
|p7|--p9| However, the stability of these steps in the pres- 
ence of hydrogen is not a priori obvious, and to our 
knowledge has not been calculated before. Therefore hy- 
drogen adsorption on three different atomic step topolo- 
gies, which have previously been discussed in context 
with th e clea. n surface, |^ will be considered in Sections 
IIIB -HID. A comparison between these results can 



lead to a better understanding of the mechanism lead- 
ing to the different chemisorption energies. Finally, the 
interaction between neighboring mono- hydri de groups on 
the surface will be discussed in Section III E . A quantita- 
tive comparison to Raschke and Hofer's [p3[ experimental 



data is given in Section III F 



II. COMPUTATIONAL METHOD 

The total-energy minimizations and geometry opti- 
mizations have been carried out using the electronic- 



structure code f hi96md l^^l in a version parallelized for 
the CRAY T3E architecture. In this code total ener- 
gies and forces acting on the atoms are calculated within 
DFT. The generalized gradient approximation (GGA) 
by Perdew and Wang H^ is applied to the exchange- 
correlation energy functional. The GGA is known to be 
distinctly more reliable than the local density approxima- 
tion (LDA) especially with respect to binding energies, 
while the LDA has proven insufhcient for the calcula- 
tion of energy barriers of reactions involving hydrogen. 
pl|-p3 1 Though recently Nachtigall and Jordan have 
argued that even the GGA is not sufficiently accurate for 
the calculation of energy barriers, the use of the GGA 
is expected to be fully adequate for the less demanding 
problem of calculating energy differences for hydrogen 
atoms adsorbed on - in our case even chemically rather 
similar - sites on the Si surface. 

The Si atoms are represented by pseudopotentials, 
which have been constructed according to Hamann's 
scheme, [ p4p5| ] consistently using the same GGA for the 
construction of the pseudopotential and the solid state 
calculations. For the hydrogen atoms we take the 
full 1 /r Coulomb potential. 

The electronic wave functions are expanded into plane 
waves up to a cut-off energy i?cut • The integration over 
the Brillouin zone is replaced by a summation over one 
or two special k points. The k points are chosen 
equidistant and aligned along the direction of the step 
edge. We assume no symmetry restrictions to the atomic 
geometry during relaxation, the only symmetry exploited 
in these calculations is time reversal. To give the reader 
an impression of the residual error due to the finite cut- 
off energy and k-point set, we will quote below numer- 
ical results at two levels of accuracy: The (30 Ry, 1 k) 
data result from geometry relaxation runs using Ecut =30 
Ry (408 eV) plane-wave cut-off energy and one special 
k point (not the F point). Important calculations have 
been repeated with 50 Ry (680 eV) cut-off energy and 2 k 
points, using the frozen geometry from the previous (30 
Ry, 1 k) run. At 50 Ry cut-off energy there are about 
153000 complex plane wave coefficients for every band 
and k-point. 

Total energies are calculated for Si(l 1 11) slabs with 
a thickness of about six atomic layers, using a super- 
cell that is periodically repeated in all three dimensions. 
This surface orientation corresponds to a miscut angle 
of 7.3° in the [110] direction away from (001). On the 
Db stepped surface the (1x2) dimerized terraces are four 
dimers wide, ensuring that we can neglect elastic step- 
step interactions when generalizing the results to 
smaller miscut angles. Furthermore, the large terrace 
width allows us to distinguish between dimers close to 
the step edge and those in the center of the terrace. The 
width of the supercell amounts to 2a, allowing for one Si- 
dimer row perpendicular to the step edge. The surface 
unit cell on the top of the slab contains either a sin- 
gle double atomic height step or a pair of single atomic 
height steps. The atoms in the bottom two Si layers 
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are fixed at their bulk positions. We use the theoretical 
equilibrium lattice constant c = a\pl = 5.450 A. On the 
bottom surface the slab is saturated with H atoms in a 
local di-hydride configuration. In case of the clean sur- 
face our supercell contains altogether 74 Si atoms and 22 
H atoms. All atoms apart from those in the bottom two 
Si layers and the bottom H-termination are allowed to re- 
lax by following the computed Hellman-Feynman forces. 
Residual forces are smaller than 3 x 10^^ eV/A. 

All results quoted below are plain total-energy differ- 
ences, not corrected for the zero-point energy (ZPE) of 
the atoms. For the vibration of the free H2 molecule the 
ZPE is fiLojl = 0.272 eV. For hydrogen atoms adsorbed 
on a Si surface in a mono-hydride like configuration the 
frequency of the Si-H vibrations does not depend much 
on the details of the configuration. Results from elec- 
tron energy loss spectroscopy indicate for Si(OOl) (2x1) 
H one stretching mode at 0.260 eV, and two bending 
modes at 0.0785 eV. |^] For two adsorbed H atoms this 
gives a ZPE of 0.417 eV. We neglect changes of the Si 
phonon frequencies due to hydrogen adsorption. Within 
this approximation the ZPE correction leads to only a 
rigid shift of all chemisorption energies by 0.145 eV. For 
deuterium this shift would be smaller by about a fac- 
tor of 1/ \/2. Differences between chemisorption energies 
are not affected by the ZPE correction, i.e., they can be 
directly read from the tables below. 



III. RESULTS AND DISCUSSION 

In the following we present our DFT results and dis- 
cuss the implications for the thermodynamics of partially 
H-covered vicinal Si surfaces. First, however, we briefly 
summarize our results for the step energies and relaxed 
geometries of the clean surface. In the subsequent Sec- 
IIID the chemisorption energies will al- 
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tions 

ways refer to the respective clean surface plus a hydrogen 
molecule at rest far away from the surface as energy zero. 
Thus, in order to find out which step topology is stable 
(at zero temperature), one first has to add the respective 
step energies from Section [II A to the chemisorption en- 



ergies before comparing the total energies of different step 
topologies with each other. As explained in the previ- 
ous Section ||, we expect the contributions from the zero 
point energy to cancel out when a total energy difference 
between rather similar, mono-hydride like configurations 
is calculated. 

It should be noted here that the probability of observ- 
ing a certain type of step (e.g., Sa-Sb or Db) along a 
step edge at some finite temperature cannot be simply 
inferred from the respective step energies per unit step 
length. The step energy difference per unit length has to 
be multiplied with an appropriate coherence length first 
(which describes the average extent of a kink-free part 
of the step) in order to obtain an excitation energy that 
makes sense to enter in a Boltzmann factor. [pOl 



Our objective is to characterize the lowest energy con- 
figurations for two H atoms adsorbed on the vicinal Si 
surface. The pairing energy for two hydrogen atoms on 
the Si(OOl) surface has been measured by Hofer, Li, and 
Heinz ||l^ to be about 0.25 eV, which is close to the 
DFT result by Northrup. This value is larger than 
the chemisorption energy difference between typical ad- 
sorption sites. Therefore we have restricted ourselves 
to configurations with two hydrogen atoms adsorbed on 
the same silicon dimer, i.e., we have only considered 
the H-paired configurations. Of course the unpaired H- 
configurations will become important when comparing to 
an experiment carried out at high temperature and low 
H-coverage, where entropy plays an essential role. We 
postpone this issue until Section IIIF. 



A. Atomic geometry of clean vicinal Si(OOl) surfaces 

A characteristic feature of the flat Si(OOl) surface are 
the Si dimers. p2| By forming these dimers the density 
of dangling bonds is halved in comparison to the bulk 
terminated surface. The Si dimers are not parallel to 
the surface. [^3|-^ In our calculation the buckling angle 
amounts to about 19° on the p(2x2) dimerized (001) ter- 
race, in agreement with experiment [ ^Jss) and previous 
DFT work. The energy lowering is driven by the 

re-hybridization of the sp^ orbitals of the Si atom that 
relaxes towards the bulk: Its bonding geometry becomes 
more planar, and therefore the bonding orbitals become 
more sp^ like, while the dangling bond gains more p char- 
acter and becomes unoccupied. The dangling bond of the 
other Si dimer atom that relaxes outwards is fully occu- 
pied. This implies that the tt bond between the Si atoms 
of the symmetric dimer is partially destroyed by buck- 
ling, however, in case of Si the energy gain due to the 
re-hybridization mechanism obviously over-compensates 
this energy cost. The buckling of the surface dimers is ac- 
companied by a change from a metallic to a semiconduct- 
ing electronic surface band structure. |^l| Due to elastical 
coupling via the atoms in the second and deeper layers 
the buckling angle alternates along the dimer row. The 
lowest energy reconstruction is p(2x2), or, even slightly 
lower, c(4x2), but the energy difference between these 
reconstructions is so small that it is irrelevant for our 
purpose. Thus we assume the p(2x2) reconstruction 
on the terraces, which allows us to keep the supercell 
small. 
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FIG. 1. Relaxed atomic geometry of single and double 
atomic height steps on the Si(OOl) surface. The orientation of 
the vicinal surface is (1 1 11), corresponding to a miscut angle 
of = 7.3°. A side view along the [110] direction parallel to 
the step edge is shown, (a) Pair of single atomic height Sa 
and Sb steps. Terraces are alternate (1x2) and (2x1) dimer- 
ized. (b) Non rebonded double atomic height Dg step, (c) 
Rebonded double atomic height Db step. In panels (a) and 
(c) the rebonded Si atoms at the Sb and Db step edge are 
denoted by arrows. 

The vicinal Si(OOl) surface can be imagined as a stair- 
case with (001) terraces separated by steps. At small 
miscut angles less than 1°, the vicinal surface con- 
sists of alternating single atomic height Sa and Sb steps, 
pO| and, consequently, the direction of dimerization ro- 
tates by 90° on successive terraces. There occurs a pre- 
sumably gradual transition initiating at about 1° to 2°, 
driven by elastic step-step interactions, and at large mis- 
cut angles the surface becomes more and more single- 
domain with predominantly double atomic height steps. 
p0| , |62 63 3^,Q These steps have been identified both ex- 
perimentally and theoretically to display the rebonded 
Db geometry, with threefold coordinated Si atoms at the 
step edge. |33j3|,|| 

We use the clean Si surface as the energy reference for 
the chemisorption energies discussed below. Therefore 
we have computed the relaxed geometries and total ener- 
gies for the three step configurations displayed in Fig. ||. 



Our results turned out to be in agreement with 
ous DFT studies of steps on the Si(OOl) surface. 
However, note that the choice of, e.g., the plane- waves 
cut-off energy is determined by the requirement to ac- 
curately describe hydrogen-containing systems, it is not 
adapted to the calculation of Si step energies, for which 
a much smaller cut-off would be sufficient. However, to 
obtain accurate step interaction energies , thicker Si slabs 
would be necessary. |48 

An Sa - Sb step pair is displayed in Fig. |l|(a). The 
(1x2) and (2x1) terraces are both 2a wide. The bonds 
between the two rebonded Si atoms and the neighbor- 
ing step-edge atom are highly strained, they are 4% and 
7% longer than the bulk Si-Si bond. The Si surface 
dimer bonds are - 2.4% shorter than the bulk bond 
length, and the buckling angle varies between 18° and 
15°, i.e., the buckling is partially suppressed close to the 
steps. The electronic structure is characterized by occu- 
pied dangling bonds at the Si dimer and step edge atoms 
that have relaxed outwards. The re-hybridization mech- 
anism works for the rebonded Si atom at the step edge 
in the same way as for the Si dimer atoms, it results in a 
height difference between neighboring rebonded step edge 
atoms of 0.7 A. This is close to the result by Boguslawski 
et al. II of 0.58 A. 

The double atomic height non-rebonded Dg and re- 
bonded Db steps are shown in Fig. |^(b) and (c), re- 
spectively. The non-rebonded geometry corresponds to a 
double layer of Si, which is terminated at the step edge. 
There are no rebonded Si atoms with highly strained 
bonds, thus the elastic interaction of the Dg steps is 
small. Among the step structures considered, this one 
has the smallest strain field on the terraces. Therefore 
the dimers in the middle of the (1x2) terrace between 
the Dg steps most closely resemble the dimers on the 
flat surface. The buckling angle of the surface dimers 
varies between 18° close to the step and 19° in the mid- 
dle of the terrace. Again, re-hybridization results in a 
pronounced buckling at the step edge. However, there 
is an important difference between the rebonded config- 
urations (a) and (c) and the non-rebonded Dg step: We 
find the Si-Si bond between the two Si atoms on the up- 
per terrace closest to the step edge to be contracted by 
3%, which indicates a strong additional 7r-bond between 
these atoms. |3| 

It was already pointed out by Chadi |^ that the non- 
rebonded Dg step is not the lowest energy configuration. 
The number of dangling bonds can be lowered by adding 
an extra "rebonded" Si atom to the D^ step edge. The 
rebonded Db structure is shown in Fig. |l|(c) , it can be de- 
scribed as a collapsed Sa - Sb step pair, with the width of 
the (2x1) terrace shrunk to zero. Similar to the Sb step, 
the bonds between the two rebonded Si atoms and the 
neighboring step-edge atoms are highly strained, 4% and 
7.7% longer than a Si-Si bulk bond (4.7% ~ 6.6% in Ref. 
pgf ) . The force dipole from these bonds induces a strain 
field and results in a considerable step-step interaction. 
We find a variation of the dimer buckling angle between 
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18° in the middle of the terrace and 16° on the lower 
terrace close to the step edge. The dimer bond lengths 
vary between 0% and 2% contraction. Compared to the 
Sa - Sb step pair we find the Db step to be ~0.02 eV 
per step-length 2a lower in energy. However, note that 
this energy difference is smaller than our rough error es- 
timate of ~0.05 eV, thus the agreement with observation 
is partially fortuitous. 

To compare the energy of the rebonded and non- 
rebonded double atomic height steps we have calculated 
the Si chemical potential by using the same supercell 
and convergence parameters as above and adding four 
Si atoms to the bulk of the slab while keeping the top 
and bottom surface structures intact. In agreement with 
previous work we find the non-rebonded Dg step to be 
energetically unfavorable. However, our energy difference 
of 0.13 eV/a is larger than Oshiyama's result of 0.06 
eV/a. This may be connected to the fact that we ob- 
serve a more pronounced relaxation of the atoms at the 
Db step edge, a height difference of the rebonded atoms 
of 0.6 A in contrast to the much smaller value of 0.17 
A in Ref. We think that our large step-edge relax- 
ation is corroborated by the similarity between the Db 
and Sb steps together with the fact that we agree with 
Boguslawski et al. with respect to the pronounced 
relaxation of the Sb step. 

In view of the smaller number of dangling bonds at 
the Db step the energy lowering with respect to the Dg 
step appears to be rather small. This can be explained 
by the extra 7r-bonds stabilizing the Dg step, and the 
large Si-Si bond strain destabilizing the Db step. Both 
effects diminish the energy difference suggested by simple 
bond-counting. 



B. Hydrogen adsorbed on surfaces with Sa-Sb steps 

The relaxed geometries for two hydrogen atoms ad- 
sorbed on a surface dimer or at the two rebonded Si 
atoms are shown in Fig. ^ for a vicinal surface with sin- 
gle atomic height steps. In panel (a) the two H atoms 
saturate the dangling bonds of the rebonded Si atom at 
the Sb step edge. Therefore there is no driving force to- 
wards buckling anymore and the two rebonded Si atoms 
become nearly equivalent. The Si-Si bond between the 
Si step edge atom and its neighbor on the upper terrace 
is strained by about 5% compared to the Si-Si bulk bond 
length. This is almost identical to the average strain of 
these bonds on the clean surface, i.e., the step interaction 
will not be much affected by adsorption. 




FIG. 2. Relaxed atomic geometries for two liydrogen atoms 
adsorbed at diflterent sites on a Si(l 1 11) surface with single 
atomic height steps. Side view along the [110] direction as in 
Fig. |l| In panels (a) - (c) the hydrogen atoms are arranged 
in a row parallel to [llO], thus only the atom in the front 
is visible, (a) Two hydrogen atoms adsorbed at the Sb step 
edge in a local mono-hydride configuration, (b) Two hydrogen 
atoms adsorbed at the Si dimer on the upper (1x2) terrace 
close to the Sb step, (c) Two hydrogen atoms adsorbed at 
the Si dimer on the lower (1x2) terrace close to the Sa step, 
(d) and (e) Two hydrogen atoms adsorbed at a dimer on the 
(2x1) terrace. 

On the whole, the geometry changes upon hydrogen 
adsorption are quite local. As can be seen in panels (b) 
- (e), the H-saturated dimers are parallel to the (001) 
surface. The respective Si-Si dimer bond lengths are 
1.6% - 2.6% larger than the Si-Si bulk bond length. The 
Si-H bond length amounts to 1.51 A. 





adsorption site 


chemisorption energy [eV] 


(a) 


H at rebonded Sb 


-2.10 (-2.13) 


(b) 


H on (1x2) terrace 


-1.78 (-1.85) 


(c) 


H on (1x2) terrace 


-1.96 (-2.01) 


(d) 


H on (2x1) terrace 


-1.90 (-1.91) 


(e) 


H on (2x1) terrace 


-1.89 (-1.89) 



TABLE I. Chemisorption energy for two H atoms ad- 
sorbed on a Si(l 1 11) surface with single atomic height steps. 
The respective adsorption geometries are shown in Fig. |^. The 
energies refer to a free H2 molecule at rest far away from an 
Sa-Sb vicinal surface, the ZPE correction is not included. 
The data in parenthesis have been calculated at 30 Ry cut-ofi^ 
energy with one special k point. The final energies have been 
calculated at 50 Ry, using two special k points in the irre- 
ducible part of the Brillouin zone, and taking the frozen ge- 
ometries from the 30 Ry, one k-point run. 
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The various chemisorption energies are summarized in 
Tab. |. We find the adsorption at the rebonded Sb step- 
edge atom to be energeticahy favored, with an energy 
difference of at least 0.14 eV with respect to adsorption 
on terrace sites. Note that for a comparison of the total 
energies of the final state configurations after hydrogen 
adsorption on, e.g., the Sa - Sb and the Db stepped 
surface, the step energy difference of 0.02eV per step- 
length 2a has to be added to the numbers given in Tab. |. 



C. Hydrogen adsorbed on surfaces with Db steps 

The relaxed geometries for two hydrogen atoms ad- 
sorbed on a vicinal silicon surface with Db steps are 
shown in Fig. ||. The observed relaxations as well as the 
chemisorption energies summarized in Tab. |l| are similar 
to the results we found for the single atomic height steps. 

In Fig. ^(a) the two hydrogen atoms saturate the dan- 
gling bonds of the two rebonded Si atoms at the step 
edge. Consequently, the buckling of the step-edge atoms 
disappears. The strain of the bond to the neighboring 
Si atom on the upper terrace amounts to 5.3% elonga- 
tion, close to the average strain of these bonds on the 
clean surface. The buckling of the surface dimers on the 
(1x2) terrace remains almost unaffected by hydrogen ad- 
sorption at the step edge. Hydrogen adsorption on the 
terrace (Fig. ^(b)-(e)) leads to local mono-hydride config- 
urations. The H-saturated Si-dimer de-buckles, and the 
Si-Si dimer bond- length expands by 1.9% - 2.5%. 



(b) H on (1x2) terrace, pos. 1 

(c) H on (1x2) terrace, pos. 2 

(d) H on (1x2) terrace, pos. 3 

(e) H on (1x2) terrace, pos. 4 



-1.77 (-1.82) 
-1.92 (-1.95) 
-1.97 (-2.00) 
-1.97 (-2.01) 



TABLE II. Chemisorption energy for two H atoms ad- 
sorbed on a Si(l 1 11) surface with double atomic height 
Db steps. The respective adsorption geometries are shown 
in Fig. ^. The energies refer to a free H2 molecule at rest 
far away from a clean DB-stepped vicinal surface, the ZPE 
correction is not included. The data in parenthesis have been 
calculated at 30 Ry cut-off energy with one special k point. 
The final energies have been calculated at 50 Ry, using two 
special k points in the irreducible part of the Brillouin zone, 
and taking the frozen geometries from the 30 Ry, one k-point 



adsorption site 


chemisorption energy [eV] 


H at rebonded Db 


-2.11 


H on (1x2) terrace, pos. 1 


-1.80 


H on (1x2) terrace, pos. 4 


-1.97 



TABLE III. Chemisorption energy for two H atoms ad- 
sorbed on a Si(117) surface with double atomic height Db 
steps. The energies refer to a free II2 molecule at rest far 
away from a clean Da-stepped Si(117) surface, the ZPE cor- 
rection is not included. The data have been calculated at 30 
Ry cut-off energy with two special k points. 




FIG. 3. Relaxed atomic geometries for two hydrogen atoms 
adsorbed at different sites on a Si(l 1 11) surface with Db 
steps. Side view along the [llO] direction as in Fig^|l[ The 
hydrogen atoms are arranged in a row parallel to [110], thus 
only the atom in the front is visible, (a) Two hydrogen atoms 
adsorbed at the Db step edge in a local mono-hydride con- 
figuration, (b) - (e) Two hydrogen atoms adsorbed at the 
various inequivalent Si dimers on the (1x2) terrace. 



adsorption site 



chemisorption energy [eV] 



(a) H at rebonded Db 



-2.09 (-2.11) 
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From the chemisorption energies in Tab. U we read 
that the H atoms preferentially adsorb on the rebonded 
Si atom at the Db step edge. In comparison to the ad- 
sorption sites on the terrace the step is energetically fa- 
vorable by at least 0.12 eV. We speculate that this may 
be due to both the different clastic relaxation energies 
at the step edge and on the terrace and some residual tt 
bond between the surface dimer atoms as opposed to the 
step edge atoms. A rigorous quantitative analysis of the 
energy contributed by the different possible mechanisms, 
however, is beyond the scope of this paper. 

Comparing the total energies for hydrogen attached 
to the Db and Sa - Sb step edge we find no significant 
energy difference within the accuracy of our approach. 

We have repeated our calculations for a smaller super- 
cell describing a Si(117) surface with Db steps, which has 
two Si dimers per (1x2) terrac e. T he chemisorption en- 
ergies are summarized in Tab. III. They compare very 
well with the values for the Si(l 1 11) surface in Tab. ||, 
i.e., we do not find any pronounced effect of miscut an- 
gle on chemisorption energies. This lends support to the 
use of the smaller, and thus computationally more con- 
venient, Si(117) supercell for the calculation of, e.g., the 
hydrogen adsorption barriers on vicinal Si(OOl) surfaces. 

A further issue that needs to be investigated is the ex- 
istence of different surface geometries corresponding to 
local minima of the total energy in configuration space. 
For this purpose we have taken the configuration dis- 
played in Fig. ^(d) and flipped the buckling angle of all 
dimers to the left hand side or to the right hand side of 
the mono-hydride group up to the neighboring step edge 
from +Q to — O, and vice versa. Then these starting con- 
figurations were relaxed in the standard way. In this way 
we found new stable configurations, which, however, are 
all energetically degenerate within the accuracy of our 
approach. At room temperature and above, the dimers 
flip rapidly between the two stable orientations. 

The large strain of the bond between the Db step-edge 
atom (i.e., the rebonded Si atom) and the neighboring 
Si atom on the upper terrace supports the speculation 
that this bond could break upon hydrogen adsorption 
and that a local di-hydride group might form at the po- 
sition of the former step edge-atom. Therefore we have 
calculated the total energy of the relaxed di-hydride ge- 
ometry on a Si(117) surface, however, we found it to be 
energetically distinctly unfavorable. The respective hy- 
drogen chemisorption energy amounts to only -0.92 eV. 



D. Hydrogen adsorbed on surfaces with 
non-rebonded Dg steps 



bonded steps discussed in the previous Sections. In par- 
ticular, the hydrogen adsorption at the step edge (i.e., on 
the adsorption sites shown in panels (a) and (b) of Fig. ^) 
is distinctly disfavored. We attribute this to the break- 
ing of the strong vr bonds between the step edge Si atom 
and the neighboring Si surface-dimer atom on the up- 
per terrace. This bond is contracted by 3% on the clean 
surface and in the configurations (c) - (f), while this con- 
traction disappears upon hydrogen adsorption on either 
of the two sites shown in Fig. ^(a) and (b) . 




FIG. 4. Relaxed atomic geometries for two hydrogen atoms 
adsorbed at different sites on a Si(l 1 11) surface with 
non-rebonded Dg steps. Side view along the [110] direction 
as in Fig. Jl] The hydrogen atoms are arranged in a row par- 
allel to [110], thus only the atom in the front is visible, (a) 
Two hydrogen atoms adsorbed at the Dg step edge in a local 
mono-hydride configuration, (b) - (f) Two hydrogen atoms 
adsorbed on the various inequivalent Si dimers on the (1x2) 
terrace. 



adsorption site 



chemisorption energy [eV] 



(a) H at non-rebonded Dg 

(b) H on (1x2) terrace, pos. 1 

(c) H on (1x2) terrace, pos. 2 

(d) H on (1x2) terrace, pos. 3 

(e) H on (1x2) terrace, pos. 4 

(f) H on (1x2) terrace, pos. 5 



-1.53 
-1.27 
-1.93 
-1.93 
-1.93 
-2.02 



TABLE IV. Chemisorption energy for two H atoms ad- 
sorbed on a Si(l 1 11) surface with non-rebonded Dg steps. 
The respective adsorption geometries are shown in Fig. ^. 
The energies refer to a free H2 molecule at rest far away from 
a vicinal Si surface with non-rebonded Dg steps. The ZPE 
correction is not included. The calculations have been carried 
out at 30 Ry, with one special k point. 



Relaxed geometries for two hydrogen atoms adsorbed 
on a Si(l 1 11) surface with non-rebonded Dg steps 
are displayed in Fig. ^, and the respective chemisorp- 
tion energies are given in Tab. fv\ . Interestingly, the 
non-rebonded step behaves totally different from the re- 
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On the whole, the absolute values of the chemisorption 
energies are smaller on the Dg stepped surface than on 
the Db stepped surface. Together with the fact that al- 
ready the clean Dg surface is energetically unfavorable 
we conclude, that the Dg step configuration is unstable 
also after hydrogen adsorption. This holds true unless 
the hydrogen coverage becomes so large that the hydro- 
gen atoms cannot be accommodated anymore on the Db 
stepped surface. Our conclusion is that after adsorption 
of two hydrogen atoms on a vicinal Si surface the most 
stable, i.e. lowest energy, configurations consist of the 
two H atoms bound to either the rebonded Sb or the 
rebonded Db step edge. 

E. Partially hydrogen-covered surfaces 

Next we consider finite hydrogen coverage and interac- 
tion effects between neighboring mono-hydride H-Si-Si- 
H groups on the surface. This will allow us to investigate 
a possible tendency towards hydrogen island formation, 
and the energetical competition of such islands, acting 
as hydrogen sinks, with the step-edge adsorption-sites 
discussed in the previous Sections. To gain further prin- 
cipal insight into the energetics of hydrogen adsorption, 
we consider both the rebonded Db and the non-rebonded 
Z?B step. In this Section we still assume the limit of large 
pairing energy, i.e., the H atoms are assumed to always 
form complete H-Si-Si-H groups. 

For the description of our approach we will restrict our- 
selves to the rebonded Db step. The different adsorption 
sites are enumerated from to 4 as in Tab. |l[ and the 
index denotes the two rebonded Si atoms at the step 
edge. The H2 chemisorption energies are taken from 
the (30 Ry cut-off energy, 1 k point) results in Tab. ||. 
This has been done for the sake of a simple compari- 
son to the results for the Dg step, for which only (30 
Ry cut-off energy, 1 k point) results are available, and 
for the sake of less expensive ab initio calculations. The 
numbers ni,i = 0, ...,4, can take the values or 1, with 
denoting a clean Si surface dimer, and 1 denoting a 
Si surface dimer with two adsorbed H atoms. We as- 
sume the following nearest-neighbor-interaction Hamil- 
tonian to describe the total chemisorption energy (take 
index 5 to be identical to index 0): 

H(■^-^=j2n^ef^^+n.n,+,w^/'-\ (1) 

i=0 

As a further approximation we assume the interaction 
parameters Wi to be equal on the terrace, i.e., wi(Db) = 
W2(Db) — u'3(Db). Effects omitted in this Hamiltonian 
are (i) the breaking- up of H-Si-Si-H pairs into single 
adsorbed H atoms, (ii) the excitation of the Db step 
into Sa - Sb step pairs, and (Hi) any effect connected 
with the configurational entropy due to dimer flips. The 
three independent interaction parameters Wi have been 



determined within DFT (at 30 Ry cut-off energy, with 
1 k point), by calculating two geometries with four ad- 
sorbed H atoms and the fully H-covered surface. The re- 
sults for the rebonded Db step are: wo(Db) = — 0.14eV, 
wi(Db) = -0.04eV, W4(Db) = -0.009eV, and sim- 
ilar calculations for the non-rebonded Dg step yield: 
wo(D^) = -0.75eV, wiiD'^) = -0.07eV, W5{r>'^) = 
— 0.007eV. Obviously, the nearest-neighbor interaction 
is attractive. We attribute it mostly to elastic interac- 
tion via the second and deeper layer Si atoms. The ex- 
traordinarily large value of |?«o(Dg)| accounts for the fact 
that two Si-Si tt bonds have already been broken by the 
first H-pair (configurations Fig. ^(a) and (b)), thereby 
favoring the adsorption of the second H-atom pair on the 
neighboring Si atoms. 

Interesting quantities that can easily be calculated 
from Eq. (|^) are the chemisorption "hole" energies, i.e., 
the chemisorption energy of the last H2 molecule that 
is finally completing the full coverage of one monolayer. 
We label the adsorption sites in analogy to Figs. ^ and 
^; to obtain the initial configuration before adsorption 
just invert the surface H-occupation (i.e., every Si dan- 
gling bond becomes H-covered, apart from those two Si 
dangling bonds that are saturated by H atoms in the re- 
spective panel of Figs. |^ or ^). The chemi sorp tion "hole" 
energies are summarized in Tabs. and |Vl| for Db and 
Dg stepped surfaces, respectively. The accuracy of our 
model Hamiltonian is demonstrated by the good agree- 
ment with two DFT test calculations, in which the atomic 
geometries have been fully relaxed. Contrary to adsorp- 
tion on the clean surfaces studied above, we find the re- 
bonded and non-rebonded step configurations to behave 
very similar in case of the "hole" energies. In particu- 
lar, for both steps the absolute value of the chemisorp- 
tion energy takes its maximum for adsorption at the step 
edge. Furthermore, the chemisorption energies on the 
terraces in Tabs. and ^ are all quite similar and be- 
low (i.e., more negative than) the values for the clean 
surfaces. This simplicity is due to two reasons: First, 
the TT bond breaking at the Dg step does not play a role 
for the hole energies. Second, there is an elastic couphng 
between the relaxation of the buckled dimers on the sur- 
face and the step edge atoms on the clean surface, which 
apparently becomes less important when the Si atoms 
are H-covered and the surface Si atoms are not buck- 
led any more. The hydrogen adsorption into the "hole" 
configuration is more exothermic than the adsorption on 
the clean surface due to the different elastic interactions: 
The dimer buckling on the clean surface tends to stabilize 
the buckled dimer and thus to de-stabilize the symmetric 
dimer of the local mono-hydride group. Our results are 
compatible with the energy difference of 0.05 eV/dimer 
between the asymmetric p(2xl) and the p(2x2) recon- 
struction of the Si(OOl) surface calculated by Ramstad et 
al. 
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adsorption site 



chemisorption energy [eV] 



(b) 
(c) 
(d) 
(e) 



rebonded Db step edge 
(1x2) terrace, pos. 1 
(1x2) terrace, pos. 2 
(1x2) terrace, pos. 3 
(1x2) terrace, pos. 4 



-2.26 (-2.25) 

-2.00 

-2.03 

-2.08 (-2.07) 
-2.06 



TABLE V. Chemisorption "hole" energy for an almost 
fuUy hydrogen-covered Si(l 1 11) surface with Db steps. The 
energies have been derived from the model Hamiltonian in 
Eq. (|l|). For comparison, the energy values in parenthesis are 
from ab initio total-energy calculations (30 Ry cut-off energy, 
1 k point). For the meaning of the labels denoting the position 
of the hole see Fig. ^. 



adsorption site 



chemisorption energy [eV] 



(a) non-rebonded Db step edge 

(b) (1x2) terrace, pos. 1 

(c) (1x2) terrace, pos. 2 

(d) (1x2) terrace, pos. 3 

(e) (1x2) terrace, pos. 4 

(f) (1x2) terrace, pos. 5 



-2.28 
-2.09 
-2.07 
-2.07 
-2.07 
-2.10 



TABLE VI. Chemisorption "hole" energy for an almost 
fully hydrogen-covered Si(l 1 11) surface with non-rebonded 
Db steps. The energies have been derived from a model 
Hamiltonian similar to Eq. (|l|). For the meaning of the la- 
bels denoting the position of the hole see Fig. ^ 



Up to this point we have focused on the chemisorption 
energies as a function of the different adsorption sites, 
and therefore it was sufhcient to consider only a single 
dimer string. In the following we will shortly discuss 
the "thermodynamic ground state" of the whole surface 
as a function of hydrogen coverage. This means that, 
while we still stick to the simple Hamiltonian in Eq. (|^) 
to describe the energy of every single dimer row on the 
Db stepped surface, we additionally allow for the ex- 
change of hydrogen atom pairs between adsorption sites 
on different dimer rows. In this way chemisorption con- 
figurations with inhomogeneous hydrogen coverage are 
included in the competition for the lowest total energy. 
In Fig. ^ the minimum of the energy given by Eq. (|l|) 
with respect to all configurations (no,rti, ...,714), rii — 
or 1, is plotted versus the number of H-Si-Si-H pairs 
n — uq + ... + n^. The full line denotes the convex hull 
and represents the lowest possible energy as a function of 
n, also accounting for the simultaneous occurrence of dif- 
ferent chemisorption configurations on the surface. Fig. ^ 
leads us to the important conclusion that the vicinal sur- 
face where all the dangling bonds of the rebonded Si 
atoms at the Db step edge are saturated by H atoms 
(configuration (no, ni, 714) = (1,0,0,0,0)) in fact rep- 
resents the ground state at this certain coverage (i.e., 
& — 1/5 for the Si(l 1 11) surface). Obviously, the in- 
teraction included in Eq. (m) is too weak to stabilize any 
other state with hydrogen "islands" on the terrace. 



^ -0.1 



c 



X 




12 3 4 

number of H-Si-Si-H groups 



9 



FIG. 5. The dots represent the minimum energy 
mm{H{no,ni, ...,714)) — fm with respect to all configurations 
(no, ni, 714) at fixed n versus the number of H-Si-Si-H 
groups n = X/i-o'^*' ^^^^ ^ surface with Db steps. 

The labels denote the lowest energy configurations. A chem- 
ical potential term /in has been subtracted for convenience, 
with fi chosen equal to the average single particle chemisorp- 
tion energy of -1.98 eV. The model Hamiltonian H from 
Eq. (|^) is used together with the ab initio energy parame- 
ters calculated with 1 k point at 30 Ry cut-off energy. The 
full line is the convex hull, which represents the thermody- 
namic ground state energy (for temperature T — > OK) as a 
function of n. 

In Fig. ^ the lowest energy configurations are shown 
for the two situations when one or two H-Si-Si-H groups 
are displaced from the Db step edge onto the terrace. In 
panel (a) the hydrogen group is forced to remain on its 
dimer string, thus the energy difference can be immedi- 
ately read from Tab. || to be equal to 0.10 eV (here we 
consistently use the results for 30 Ry, 1 k point; for 50 
Ry, 2 k points the energy difference would amount to 0.12 
eV). The energy difference becomes slightly smaller, 0.09 
eV, for the configuration displayed in panel (b), when 
the hydrogen atoms are transferred onto another dimer 
string. Finally, when two hydrogen groups arc displaced 
as shown in panel (c), the energy difference amounts to 
0.16 eV. 





(a) I I I I I (b) \\ \]\ (c) 

FIG. 6. Metastable configurations for hydrogen adsorbed 
on a Si(l 1 11) surface with Db steps and hydrogen cov- 
erage 1/5. The Si surface dimers and the step-edge atoms 
are denoted by lines. On the terrace the dimers arrange in 
dimer strings. The open and filled circles denote the position 
of the mono-hydride H-Si-Si-H groups in the ground state 
and in the excited state configuration, respectively, (a) One 
mono-hydride group displaced from the step edge onto the 
terrace within one dimer row. (b) One mono-hydride group 
displaced from the step edge onto a terrace site of another 
dimer row. (c) Two mono-hydride groups displaced onto an- 
other single dimer row. 



F. Thermodynamics at low hydrogen coverage 

In this Section we will present a comparison between 
our ab initio results and the experimental hydrogen ad- 
sorption data by Raschke and Hofer described in the pre- 
ceding paper. |Q They have saturated the Db step edges 
of vicinal Si(OOl) surfaces miscut by 2.5° and 5.5° to- 
wards the [110] direction with molecular hydrogen at a 
surface temperature, which is sufficiently low to suppress 
the diffusion of the hydrogen atoms onto the terrace. 



This non-equilibrium configuration can be prepared, be- 
cause there is a considerable energy barrier towards the 
dissociative adsorption of the II2 molecules onto the fiat 
surface or the terrace, while the adsorption-energy bar- 
rier at the step edge is very small. |3^,^ At temper- 
atures between 618 K - 680 K partial thermal equilib- 
rium between the adsorption sites at the step and on the 
terrace is established. After equilibration the configu- 
ration is frozen-in by rapidly cooling down the sample, 
and the residual hydrogen coverage at the Db step-edge 
site is measured. This is accomplished by determining 
the amount of hydrogen, expressed as a change of cov- 
erage A0, necessary to saturate the depleted Si step- 
edge dangling-bonds again. The resulting step-edge H- 
occupation-probability po = 1 — ^©/©itopi based on the 
experimental step-edge H-saturation coverage 6|tip, is 
denoted by the filled circles and the square in Figs, j^and 




200 400 600 

temperature [K] 



800 



FIG. 7. Probability po, that the Si dangling bond at the 
Db step edge is saturated with a hydrogen atom, versus the 
surface temperature for two vicinal surfaces with miscut an- 
gles of 2.45° and 5.4° degrees. The total hydrogen coverage 
amounts to two hydrogen atoms per dimer string, i.e., it has 
been chosen equal to 1/{N + 1) independent of temperature. 
Different line styles are used to distinguish among the various 
approximations. Dotted line: Hamiltonian Eq. (|l[), including 
the nearest neighbor interaction term. Dashed line: Same 
Hamiltonian, but all interactions Wi set to zero. Full line: 
Hamiltonian Eq. (^) with pairing energy tpair = 0.25eV for 
all sites. The symbols denote experimental data by Raschke 
and Hofer for vicinal surfaces with a miscut angle of 2.5° (cir- 
cles) and 5.5° (square). 

To compare to experiment we generalize the Hamilto- 
nian in Eq. (|l|) from the Si(l 1 11) surface with N—A Si 
dimers per terrace to DB-stepped vicinal surfaces with 
smaller miscut angle and iV>4 dimers per terrace. This 
is done by simply replicating the energy corresponding to 
configuration (c) on the Db stepped surface (see Fig. ||) 
together with the interaction parameter ^^(Db), which is 
assumed to be constant on the terrace. The configuration 
(c) was chosen, because its chemisorption energy is clos- 
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est to the chemisorption energy on the terraces bounded 
by Dg steps. As aheady noted above, the strain field of 
the non-rebonded Dg step is distinctly weaker than that 
of the rebonded Db step. Therefore the terraces on the 
Dg stepped surface resemble the flat surface (or, equiv- 
alently, the middle part of a wide terrace) more closely 
than any of the other configurations. For the miscut an- 
gles 5.4° and 2.45°, corresponding to N=6 and A^=15 
dimers per terrace, we have calculated the step-edge oc- 
cupation probability po = ('^o) at fixed hydrogen cover- 
age Q — 1/{N + 1) as a function of surface temperature 
within the grand canonical ensemble, see the dotted line 
in Fig. !?[ At this small hydrogen coverage the interaction 
term in the Hamiltonian (|l]) does not play any important 
role, as can be judged by comparison to the dashed line 
in Fig. 1^, which was calculated in exactly the same way, 
but with all interaction parameters Wi set to zero. 

However, it is well known that at small hydro- 
gen coverage configurational entropy dominates over the 
hydrogen pairing energy. Single, un-paired chemisorbed 
H-atoms become more and more frequent at low EE- 
coverage. To account for this process we generalize our 
model Hamiltonian. The occupation numbers nf and 
nf can take the values or 1, depending on whether the 
Si atom "A" or "_B" of the i-th Si surface-dimer is H- 
saturated or not. For simplicity we omit the interaction 
between H-atoms on different dimers, i.e., we set all Wi 
to zero. Furthermore we assume a site-independent pair- 
ing energy epair of 0.25 eV, in agreement with previous 
experimental and theoretical work for the flat surface. 
The Hamiltonian is: 

i=0 

+ ^(nf(l-"f)+"f(l-"f))W- (2) 

We have calculated the H-occupation probability at the 
Db step edge, po = {uq + Uq) /2 for this model. The re- 
sults are plotted as full lines in Fig. 0. Among the model 
Hamiltonians discussed above we expect this approach 
to yield the most realistic description of the Si surface 
in thermodynamic equilibrium at low H-coverage. The- 
ory and experiment agree with respect to the fact that 
the hydrogen more tightly binds to the step edge than 
on the terrace. However, as can be read from Fig. 
there remains a slight quantitative discrepancy. Because 
the chemisorption energies at the Db and at the Sb step 
edge are similar, we do not attribute this difference to 
the existence of single atomic height steps on the ex- 
perimental sample. Instead, we argue that at least part 
of the difference may be due to the incomplete thermal 
equilibration of the surface in experiment. There is quali- 
tative agreement between theory |65|-|67| and experiment 
p8[ that the diffusion of hydrogen atoms on the Si(OOl) 
surface is highly anisotropic and fast along the dimer 
rows. The experimental activation energy for hop- 
ping along a dimer row amounts to Eb = 1-Q8 eV when 



an attempt frequency of i'o=10^'^ s~^ is assumed. At 
a surface temperature of T=620 K this corresponds to a 
hopping rate ly = vq exp{— Eb / ksT) —0.2 s^^, i.e., many 
of these diffusion events occur during the time interval of 
the order of lO'^ s during which the surface is kept at el- 
evated temperature in Hofer and Raschke's experiment. 
Similarly, the intra-dimer hopping is predicted to be as 
fast as the intra-row hopping of the hydrogen atoms. ]6^ ] 
Thus we may safely assume a perfect thermal scrambling 
of the H-atom adsorption sites within a single dimer row. 
The hydrogen inter-row diffusion, on the other hand, is 
distinctly slower. The theoretical barriers vary from 1.8 
eV Q , which translates into a hopping rate of roughly 
0.02 s~^, to barriers much larger than 2 eV |^,|^, where 
hopping perpendicular to the dimer rows would be almost 
completely suppressed. To investigate this issue we have 
calculated the thermal expectation value {n^ + n^)/2 for 
the hydrogen coverage at the step edge from the Hamil- 
tonian (^) within the canonical ensemble. In this case the 
number of hydrogen atoms on every single dimer string 
is restricted to be equal to 2. The result in Fig. [s] shows 
a better agreement with Raschke and Hofer's experiment 
than Fig. Ui i.e., our calculations are consistent with the 
notion that at 620 K - 680 K the hopping of H atoms 
across the row and the exchange of H atoms between 
rows on different terraces is indeed a rare process within 
a typical experimental time interval of 10^ - 10'^ s. 




temperature [K] 

FIG. 8. Same as Fig. ^. However, in this case only partial 
thermal equilibration of the system is assumed: The adsorbed 
hydrogen atoms are allowed to diffuse only along the dimer 
row, or to hop from one Si atom to the other Si atom within 
the same dimer. Results are based on the Hamiltonian in 
Eq. ^ and a canonical ensemble. The symbols denote ex- 
perimental data by Raschke and Hofer for vicinal surfaces 
with a miscut angle of 2.5° (circles) and 5.5° (square). 



IV. SUMMARY AND CONCLUSIONS 

The chemisorption energies for H2 on vicinal Si(OOl) 
surfaces have been calculated using a well established ab 
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initio total-energy technique. |4C| ] The various adsorp- 
tion sites on the terraces and at the step edges have been 
compared to each other. Both the single atomic height 
Sa-Sb steps and the rebonded and non-rebonded dou- 
ble atomic height Db and Dg steps have been consid- 
ered. The results confirm recent experimental findings 
by Raschke and Hofer for the adsorption of molecu- 
lar hydrogen on vicinal Si surfaces. 

In particular we conclude from our calculations: 
(i) On a surface with Db steps hydrogen preferentially 
binds to the step-edge Si atoms. The chemisorption en- 
ergy difference with respect to the terrace sites is > 0.12 
eV per H2 molecule. 

{ii) Adsorption on a vicinal surface with single atomic 
height steps is similar to adsorption on the Db stepped 
surface. The hydrogen atoms preferentially bind to the 
Sb step edge, with nearly the same chemisorption energy 
as for the Db step edge. 

(Hi) The non-rebonded Dg step, which is known to be 
unstable on the clean surface, [ p7p9| ] remains unstable 
also after hydrogen adsorption. 

(iv) The interaction between neighboring mono-hydride 
groups on the Si surface has been calculated and was 
found to be locally attractive. However, this interaction 
is too weak to be of importance for the conclusions drawn 
in this paper. 

(w) The hydrogen coverage of the Db step-edge has been 
calculated as function of temperature. We assume that 
the diffusing hydrogen atoms are confined to a single 
dimer row and obtain qualitative and even quantitative 
agreement with the experiment. |Q 

This work supplies the basis for a theoretical investi- 
gation of the H2 adsorption energy barriers on vicinal 
Si(OOl) surfaces, especially for a comparison between the 
dissociative adsorption at the step edge and on the Si 
dimers of the terrace. 
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